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Abstract: 

We argue that parton distributions in coordinate space provide a more natural object for 
nonperturbative methods compared to the usual momentum distributions in which the 
physics of different longitudinal distances is being mixed. To illustrate the advantages 
of the coordinate space formulation, we calculate the coordinate space distributions for 
valence quarks in the proton using the QCD sum rule approach. A remarkable agreement 
is found between the calculated and the experimentally measured u-quark distribution 
up to light-cone distances A~ = A — A 3 of order ~ 1 fm in the proton rest frame. The 
calculation for valence d quarks gives much worse results; the reasons for this discrepancy 
are discussed. 
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1 Introduction 



Deep inelastic lepton-hadron scattering has proved to be the best testing ground for per- 
turbative QCD. Thanks to the celebrated factorization theorems KJ, which can be derived 
quite rigorously in this case by using the Operator Product Expansion (OPE), the entire 
Q 2 dependence of the cross section can be calculated perturbatively, while all dynamical 
effects of large distances are included in a set of one-particle parton distribution functions 
given at a certain reference scale. Determination of the set of partonic distributions - 
quark, antiquark, and gluon — is an ultimate goal for the experimental studies of the deep 
inelastic scattering, and also provides a challenging task for nonperturbative approaches 
to QCD. 

In the past ten years remarkable progress has been made at experimental side, and 
apart from the region of small Bjorken x, there is not much controversy regarding the 
existing parametrizations of parton distributions. The theoretical progress has been much 
more moderate. Apart from several quark-model or MIT bag model calculations, there 
have been relatively few attempts to determine parton distributions e.g., from QCD sum 
rules. The problem has proved to be difficult for the theory. Purpose of this paper 
is to point out that major part of theoretical problems in the calculations of parton 
distributions is due to the fact that distribution functions in momentum space for each 
particular value of the momentum fraction x receive contributions from both small and 
large longitudinal distances, which correspond to different physics and are difficult to treat 
simultaneously. We argue that longitudinal distance distributions are much easier from 
the theoretical point of view than momentum distributions. They can be extracted from 
data with marginal complications. In what follows we shall demonstrate that working 
with coordinate space distributions gives us a selfconsistent formalism, which is not more 
complicated than the standard one and its relation to the OPE is much closer. We 
illlustrate the advantages of this approach by the calculation of valence quark distributions 
in the framework of QCD sum rules. 

Obviously, the longitudinal distance distributions are simply Fourier transforms of the 
momentum distributions. For example for the valence quark distribution one gets 

QvaiO, At 2 ) = / du cos (uz)q v (u,fi 2 ), (1) 
Jo 

where u is the momentum fraction. The physical interpretation of the variable z has 
been discussed in the literature since a long time || f|. In the center-of-mass (CM) 
system of the target a deep-inelastic probe 7* (photon or neutrino) is converted into a 
quark- antiquark pair at some space-time location. At large Q 2 the pair travels with the 
speed of light along a light-like path, interacts with the nucleon and is converted back 
into the probe. The time interval between the conversion points 7* — > qq and qq —>■ 7* in 
the CM system, so called Ioffe time 77, measures the light-like distance essential for the 
process. The Lorentz-invariant variable related to 77 is denoted by z. In the CM system 
the relation between these two variables takes a simple form 77 = -h, where M is the mass 
of the target (nucleon). For simplicity we shall call z just the Ioffe time. 
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The idea to study z-distributions in order to understand the relative importance of 
various light-cone distances for deep inelastic scattering was proposed as early as in 1970 
IJ. There, the first analysis based on experimental data available at that time was 
presented. The experimental information about the Ioffe-time distributions available to- 
day is summarised in Fig. [I]. We plot Fourier transforms for two representative sets of 
parametrizations || - the NLO Gliick, Reya, Vogt parametrization || and the CTEQ 
parametrization || - for the valence quark (u and d), quark plus antiquark, and gluon 
distributions at Q 2 = 4 GeV 2 . Exact definitions are given in Sect. 2, where the connection 
of the Ioffe-time distributions and the OPE is discussed in detail. We note that all existing 
parametrizations are in a fairly good agreement at small z, while at large z there is some 
discrepancy. In order to get some insight about the transition from the momentum to the 
coordinate space, we plot in Fig. [|the Ioffe-time valence quark distribution corresponding 
to a simple ansatz 

q v ( u ) = Nu a (l-uf 

with a = —0.5 and (3 = 3, as suggested by the Regge theory and the perturbative QCD ar- 
guments, respectively. The normalisation factor iV ensures that Q va i(0) = Jq duqy{u) = 1. 
The shape of Q VSL \ depends very weakly on the exact value of the parameter (3 as long as 
a <C (3. The large-z behavior can readily be obtained from the standard expressions for 
the asymptotic expansion of hypergeometric functions.^ 

Q{z) ~ - sin (-oO-L-^ + /^cos (-a) \ a+2 +■■■ (2) 

Aa anticipated, the large-z asymptotics of Q(z) ~ z a+1 is fixed by the Regge behaviour. 
The valence quark distributions decrease at large z as ~ 1/ yfz. The sea quark and gluon 
distributions should approach a constant, or even may be rising functions of z at large z. 

An interesting question is at which values of z, or equivalently at which longitudi- 
nal distances, the behavior of parton distributions is already determined by the Regge 
asymptotics. The dashed line on Fig. [2] shows the asymptotic expansion of the function 
Qvai{ z ) which matches almost perfectly the true behaviour when z > 6. Because in the 
target (lab) frame z = \M8, where 5 denotes the distance along the light-cone, in the 
case of the nucleon these values of z correspond to the light-cone distances 5 > 2.5 fm.0 
A nonperturbative calculation of the parton distribution in this region is an extremely 
intricate theoretical problem, which is essentially equivalent to providing a nonperturba- 
tive input to the Lipatov's pomeron || (or Reggeon, in case of valence distributions). On 
the other hand, calculation of the distributions at sufficiently small z may be within reach 
of existing models, the QCD sum rules or the lattice calculations. Note that all parton 
distributions represented in Fig. [I] are very smooth at small z, which suggests that one 

t The asymptotical expansion contains in addition oscillating terms ~ l/z l3+1 but they can be ne- 
glected since a <C /3 in realistic case. 

avoid confusion, note that we are speaking here about light-cone distances, which in our definition 
are factor two larger than physical longitudinal distances. 



2 



should be able to reproduce them in this region with only few terms in the Taylor expan- 
sion around z = 0. These terms are related to the first few moments of the momentum 
distributions, in other words to nucleon expectation values of a few local operators of low 
dimension (see below). The most important question is whether there exists a "matching 
window" , where both the Regge asymptotic formulae and the small z expansion are ap- 
plicable. Provided the answer is positive, one could hope to get a quantitative description 
of the parton distributions in the whole z range, matching these two different inputs at a 
certain intermediate value of z - note similarity with the usual QCD sum rule program. 
Thus, the problem of calculating the parton distributions can be posed as a problem of 
calculating the distributions at distances of order 2 — 3 fm at the light-cone. As we shall 
see below, the standard QCD sum rules are sufficient for this purpose for the valence 
u-quark distributions, but fail for the d-quarks. 

Playing around with typical parametrizations for parton distributions which are used 
in modern experimental analysis, one can convince oneself that in all cases the onset of 
the Regge behavior corresponds to values of z ~ 5 — 8, see Fig. [|. Another useful example 
is given by the polarized gluon distribution. In this case we define 

AG(z,/i 2 )= / duusm(uz)Ag(u, /i 2 ) (3) 
Jo 

where Ag(u,fi 2 ) is the usual polarized gluon distribution depending on the momentum 
fraction u (see Sect. 2 for details). Note that the gluon polarisation 

rl roc 

Ag = / duAg{u) = / dzAG{z) (4) 



A typical shape of AG(z) is shown in Fig. |3|. It has been obtained with a simple model 
of Ag(u) = N G u a {l — u)^, and the solid curve corresponds to a = and (3 — 4. The 
normalization constant Nq is chosen in such a way that the gluon polarisation Ag = 0.5. 
The two short-dashed curves were obtained by taking (3 equal to 3.5 and 4.5 respectively, 
and keeping a = and Ag = 0.5 fixed. The corresponding variation of AG(z) is rather 
mild, and one can conclude that the behaviour of Ag(u) at small u combined with the 
value of Ag determine to large extent the shape of AG(z). By the same argument, 
knowledge of AG(z) up to the point of maximum, which is again at z ~ 6, is enough to 
estimate the value of Ag within, say, 50% accuracy. Note however that because of the 
more complicated shape of AG(z), asymptotic expansion shown by the long-dashed curve 
in Fig. ^| starts to be valid at larger values of z ~ 10. 

Our presentation is organized as follows. In Sect. 2 we discuss the theoretical frame- 
work for the introduction of parton distributions in coordinate space as matrix elements 
of nonlocal operators, and we describe their Q 2 evolution. The presentation in this section 
mainly follows Refs.|| |7]]. In Sec. 3 we give the QCD sum rule calculation of valence quark 
distributions in coordinate space, and compare our approach to the direct calculation in 



the momentum space in Ref. fl5|. Sect. 4 is reserved for a summary and conclusions. 



Some technical details of the sum rule calculation are presented in the Appendix. 
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2 Covariant Definition of Ioffe Time Distributions 
and their Q 2 Evolution 

An intuitive discussion of the space-time picture of deep inelastic scattering in late 
60-s can be put on a rigorous footing using the formalism of the Operator Product 
Expansion (OPE). We are going to demonstrate that the Ioffe time distributions arise 
naturally in this framework as reduced matrix elements of nonlocal string operators on 
the light-cone. Our presentation essentially follows Refs. [§, 0]. 

It is well known that the deep-inelastic ep scattering cross-section is related to the 
matrix element: 

T, v = J d 4 y ex V (tqy)(P | T\j^j v {Q)\ I P) (5) 

where | P) represents a proton with momentum P and j is an electromagnetic current 
operator. This quantity describes hadronic part of the process. The Operator Product 
Expansion applied to T[j(y)j(0)] gives rise to its systematic expansion in powers of the 
small parameter A 2 /Q 2 where A is the QCD scale of the order of 200 MeV and Q 2 = - q 2 
is the virtuality of the deep-inelastic probe. To the leading, twist-2 accuracy, i.e. when 
all powers of A 2 /Q 2 are neglected, the quark and gluon operators appearing in the OPE 
of T[j(y)j(0)] have the form: 

0W-Mn( O ) = If (0){7 Ml ^ 2 ...^"} gT ^(0) 

0W-^(o) = ^{G^^O)^ 2 ...^"- 1 ^^^)^ , (6) 
y 2 

where = <9 M — igA^ 1 denotes covariant derivative, and \1/ and G^ v are the quark field 
and the gluon field strength, respectively. The subscript ST denotes the symmetric and 
traceless part of a Lorentz tensor. 

The operators defined in (|J) form an irreducible representation of the Lorentz group. 
Their reduced matrix elements ((O™)) and ((0™)) are defined by the relations: 

( P | 0Mi-Mn( O ) | p) = ((O™)) {P^ . . . P^} ST 

( P | ^i-^(o) | p) = ((0™)) {P^ . . . P^} ST . (7) 

According to the standard analysis, the matrix elements are related to the moments of 
familiar quark and gluon distribution functions q(x), q(x) and g(x). 

«Oy»^ = / 1 rfxx"- 1 (g(x,/i 2 ) + (-irg(x,/i 2 )) (8) 

J 

((O n g ))^ = J^dxx^gix^ 2 ). (9) 

Equation (H) holds only for even values of n. Note that through renormalization the 
operators (j^) acquire a scale dependence, which is related to the scale dependence of the 
parton distribution in (||) and (|). 
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An alternative representation for relations (0) and (BI) has been noticed in || long ago, 
and we want to introduce it now. For that purpose let us define two light-like vectors 
such that n 2 = A 2 = n ■ A = 0. Our convention is such that and for any 
a + a 3 . The vector A is just proportional to n, A = ^Sn, where 
A 3 is the distance along the light cone. Following Ref. H] we can write 



n M and A M 
vector a, n ■ a = a + 
5 = A~ = A - A3 ; 



(P|f(A)[A;0]*(0)|P) M2 



for quark distributions and 



2(P 



du 



(q(u, /z 2 ) exp (iuz) — q(u, u 2 ) exp (— iuz) 



(10) 



(P|GUA)[A;0]G*(0)|P>^ 



for gluons. In the above formulae z 
for the path-ordered exponential: 



u = A(P-nf dug(u,u 2 )u cos (uz) (11) 
J o 

P ■ A and we have introduced the notation [A; 0] 



[A;0] = Pexp 



igA $ £ dsA^(As) 



'121 



which is necessary for gauge independence of the parton distributions considered. An easy 
way to obtain the relations fllOD, (|TTD and (p~3|) is to insert the complete set of intermediate 
light-cone quark or gluon states between the field operators at the RHS working in the 
Schwinger gauge: A-A(A) = 0. Taylor expansion in A of both sides of (0) and (0) gives 
exactly the set of relations (|8]) and (|9|) between the matrix elements of local operators 
and the moments of structure functions. Note that because n and A are light-like and 
proportional, the local operators arising here automatically are of twist 2. 

Fourier transformation of (|T0|) and ([O]) gives a gauge invariant definition of parton 
momentum distributions in terms of reduced matrix elements of leading twist nonlocal 
operators at the light-cone 0. On the other hand, it is possible to demonstrate that 
the usual program of the OPE can be formulated directly in terms of nonlocal light- 
cone operators H. Thus, this formalism is consistent. For completeness, we quote the 
definition for the polarised gluon distribution []TTJ, [12]], see discussion after Eq.@: 



(P,S|G^(A)[A;0]G*(0)|P,%. ? 



where S is the nucleon spin vector normalized by S 2 = —M 



Ai(P-n)(S-n) / du Ag(u, [i 2 )u sin (uz) , (13) 



N ■ 



Taking C-odd and C-even combinations of the LHS of (0) we arrive at the definitions 
involving C-odd (valence) and C-even combinations of parton densities: 



(P | *(A) j< [A; 0]*(0) | P)„2 + (A -> -A) = 4(P • n) \ duq v (u, u 2 ) cos (zu) (14) 



and 



(P | tf(A)^[A;0]tf(0) | P) fi 2-(A -> -A) = 4z(P-n) 



du 



q(u, u 2 ) + q(u, u 2 ) sin (zu) 



(15) 
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where qv(u,fj?) = q(u, /i 2 ) — q(u,fi 2 ), and fl = n^. 

As mentioned above, the conventional procedure is to Fourier transform the above 
formulae ending up with the parton distributions in momentum space. The main thrust 
of our paper is to point out that the matrix elements appearing on the LHS of the 
equations (|ITj|) — (P~5|) have a clear physical interpretation as the parton distributions in the 



longitudinal coordinates, and are more adequate for the application of nonperturbative 
methods, retaining at the same time the whole physical content of the momentum space 
description. We define Ioffe time distributions by: 



(P | #(A) [A; 0]*(0) | P},2 + (A - -A) = 4(P • n) Q^(z, /i 2 ) , (16) 

(P | (A) i [A; 0M0) | P), 2 - (A - - A) = U{P ■ n) Q(z, /i 2 ) (17) 
for quarks, and 

(P | G*(A)[A;0]G*(0) | P)pn»n v = 4(P ■ n) 2 G(z, fi 2 ) , (18) 

(P, S | G^(A)[A; 0}Gl(0) | P, S)^n v = 4z(P • n)(S ■ n)AG(z, fi 2 ) , (19) 
for nonpolarized and polarized gluon distributions, respectively. Comparing to (|lO|)-(p^) 



we arrive at the relations between the momentum and coordinate space distributions in 

According to the standard discussion the scale dependence of the longitudinal momen- 
tum parton distributions is governed by GLAP equations written in momentum space. 
However, it is also possible to derive the corresponding evolution directly in coordinate 
space || 0. The coordinate-space version of the LLA evolution equations has been ob- 
tained in |0] in the form of equations describing the normalization-point dependence of 
the non-local operators fllCf ) and (pTT|) . Taking the forward nucleon matrix element and 
making use of relations (JTTJ) , (|i~4|) and (|i~5|) one can derive evolution equations for Ioffe 
time distributions Q Y& \{z, fi 2 ), Q(z,fi 2 ) and G(z,n 2 ). 

To one loop accuracy the scale dependence of the valence quark Ioffe time distribution 
(16) is governed by 

(y<s Li 2 

Qvai(z;i4) = Qvai(^;/i 2 ) - —C F \n^ / duK(u)Q val (uz; /i 2 ). (20) 

2-n ii{ Jo 

The kernel K(u) is given by 

1 r?/l 

K(u) = -8(u)-u-2 - (21) 

2 lui + 

where u — 1 — u and for any function f(u) 



du 



o 



f(u)= I du-(f(u)-f(l)) (22) 
Jo u 
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In the flavour-singlet channel the evolution mixes, as expected, quark and gluon dis- 
tributions Q(z,/j, 2 ) AHD and G(^,/i 2 ) ©: 



Q(z;i4) 



where 



Q(z; n\ 
zG(z; \i\ 



■ — In ^| / du 

2n uf Jo 



C F K QQ {u) N f K QG (u) 
C F K GQ (u) N c K GG (u) 



Q(uz; u\) 
zG{uz\ uf) 



(23) 



k qq(u) 
K GG {u) 



u 



1 1 N f . , 



+ 2(u 6 - u ) 



K QG (U) 

K gq (u) = -«?(«)- 2iZ 



— m(2m 2 + 3m) 
3 



(24) 



where N G is the number of colours and Nf is the number of active flavours. Equations 
( pop and ([23|) allow for a systematic study of the evolution of Ioffe time parton densities 
exactly in the same manner as the conventional GLAP equations do for the longitudi- 
nal momentum parton densities. Indeed, let us rewrite (|20| ) and (|23"f ) as the evolution 
equations for flavour non-singlet and flavour singlet distributions: 



(9 

2 /~\ l 1 



2ix 



-C f 



duK(u)Q val (uz; /j 2 ). 



(25) 



and 



_d_ 



Q{z-^) 
zG(z; u 2 ) 



OsljJ 2 

2n 



du 



C F K Q 
CpK, 



gq{u) 



N f K QG {u) 
N G K GG (u) 



Q(uz; u 2 ) 
zG(uz; u 2 ) 



(26) 

which can be identified as the RGE for twist- 2 Ioffe time distributions. 

A beautiful feature of the RGE equations (|25| ) and (|26| ) is that they explicitely show 
the relevance of the short distance expansion: To calculate the QCD evolution of the 
distributions, one needs to know them at a certain reference scale at smaller values 
of the Ioffe time. Stated differently, the evolution equation for the nonlocal operators 
flIUD-([T3]) involves these operators at quark- ant iquark distances smaller than the initial 
separation. This is in contrast to the evolution of fragmentation functions, which is 
essentially nonlocal in the coordinate space |TtJ . 



As it is well known, the integrodifferential equations (25 ) and (|26|) can be transformed 
into ordinary differential equations by Mellin transformation 



F(z) 



F(u) 



dzz v ~ 1 F{z) 



As a consequence we obtain in the flavour non-singlet case: 



2 9 A i 2n a s(v 2 

/^Qva,(^) = - — 



-71 



l)Qvai(^;M 2 



(27) 
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where we have introduced the function 



7(1/) = 2C F [ duK(u)u v - x {2i 



Equation (|27| ) has the well known solution 



/ I 2\ \ T(-"+l))/6 

Qvai(^^)=Mf|( Qvai(^;^), (29) 

where in QCD b = ^ — ^Nf. The Ioffe time distribution at the scale /i 2 , can be obtained 
with the help of the inverse Mellin transformation: 

1 rc+ioo 

F{v) -> = — ; / dv z~ v F{y) . 

2wi Jc-ioo 

Explicitely, one gets 

y(l/2-iv)/b 



Q v&l {z- fi) = / ^ ^|( / rf M ^- 1 / 2 g va i(^; Mi) • (30) 

Typical results of the low-scale evolution i.e., in the low /i 2 range are illustrated on 
Figures [| and |5|. Figure ^ shows valence and gluon GRV Ioffe time distributions evolved 
between fi 2 = 4 and 20 GeV 2 . Figure |5] shows the same for u- and d-quark distributions. 
Note that in this range of scales the evolution affects mainly the large-z behavior of the 
distributions. 



3 Valence Quark Distributions from QCD Sum Rules 

We can now summarize our discussion in the following way. We have analysed the reduced 
matrix elements of QCD string operators of twist 2 as a function of the light-like sepa- 
ration between fields. We have found very smooth behaviour which makes such objects 
convenient for theoretical studies. It can be also shown that the large separations are 
dominated by the asymptotics of the corresponding structure function at small values of 
Bjorken x. Once this asymptotics is known from e.g. the Regge arguments, the remaining 
non-trivial information is contained in the domain of moderately large separations. This 
region could be accesible to presently developed analytical methods like QCD sum rules, 
instanton models of QCD vacuum [p3[| , or lattice calculations. 

In this paper a QCD sum rule calculation is carried out for the valence quark distri- 
butions. In the last decade the QCD sum rule approach has been applied succesfully to 
a variety of problems, including estimation of hadron masses and couplings, elastic and 
transition form factors etc. In the context of this paper it is necessary to mention the 



calculation of the fraction of proton momentum carried by gluons in Ref. |TB|, III] and the 



calculation of structure functions at intermediate values of Bjorken variable in Ref. [15 
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In what follows, we shall often refer to this latter analysis to compare the calculations in 
coordinate and momentum space. 

The basic idea of the QCD sum rule technique is to use duality between hadronic 
and partonic representations of a suitable correlation function to extract the quantity of 
interest by requiring that the two descriptions match each other at intermediate scales. 
To calculate the valence quark distributions, we choose to work with 

Tl u > d = i 2 j d 4 xd 4 yex V (ip ■ x + iq ■ y)(Q | T[rj(x)fj(0)d u /(y + ^; y - |)] | 0) (31) 
where 

o l s(y + y ; y - f ) = v(v + yhXfo + ^v~ jMv - y) + (A - -A) . (32) 

Here i = u,d denotes quark flavor, = (1, 0, 0, —1) is the "unit" light-like vector, n 2 = 0, 
and the splitting A M is light-like and proportional to n M . 

Finally, rj(x) is the standard interpolating current for the proton |16 

v (x) = e^u^xfC^u^x^d^x) (33) 

where u(x) and d(x) denote u- and d-quark fields, respectively and a, b, c are color indices. 

According to Eq. fll6| ) the proton matrix element of O l (y + ^; y — -j) defines the 
valence quark distribution. In the following we choose a special kinematics taking the 
momentum transfer to be light-like q 2 = and orthogonal to the interquark separation 
i.e., A • q — 0. In this case, the nucleon contribution to the correlation function (pi]) can 
be extracted in the form 

1 4 A 2 z" 1 

-Tr^f n"' d = —t. ^ . . N T-prrip ■ n) 2 / du qy (u) cos (up ■ A) + continuum (34) 

4 (p 2 -M 2 N )(p 2 -M 2 N y F io HV y y 

where p\ = p 2 and pi = (p + q) 2 - The coupling X N is defined by 

(n\ v (o)\p,N) = x N u N (P) 

where un{P) is the nucleon spinor. 

Note that the operator O in (|3^) is essentially the point-splitted vector current. In 
the limit A = a Ward identity relates the three-point correlation function fl3"T| ) to the 
derivative of the two-point correlation function of two nucleon currents 

IT< d (A = 0) = N u ' v f 1 dvn,^-U i2) (p + vq) , 

Jo dpf, 

n (2) (p) = i f d 4 xexp (ip ■ x)(tt | T [77(20*7(0)] I Q) , (35) 



where N u = 2, N d = 1 are the numbers of valence quarks in the proton. The derivation 
of fl3"5]) and of a more general Ward identity for arbitrary separation A is given in the 
Appendix. Note that 

-Tr^II^Qo) = 9 — — -(p ■ n) + continuum. (36) 
4 ' Mfj — p 1 
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Substituting ( |36| ) in ( |35|) and comparing to (B3I) one obtains the normalization conditions 



C/vai(0) = [ 1 duq u (u)=N u = 2, 
Jo 

^vai(O) = [ 1 duq d (u) = N d = l, (37) 



which are exact in the QCD sum rule approach, provided the Ward identity is not spoiled 
by the continuum subtraction (see below). 

Main task is the calculation of the correlation function (|3l|) in QCD. If both p 2 and q 2 
are sufficiently large (and negative) the dominant contributions come from small distances 
x and y of order l/—p 2 and l/—q 2 , respectively. Thus the standard machinery of the short- 
distance expansion is applicable, allowing to express the result as a power series in terms 
of vacuum quark and gluon condensates. In the case of the forward matrix elements (i.e., 
for q 2 = 0) the situation is more involved because the relevant distances in the t-channel 
are not constrained by the external momenta and can be arbitrarily large. 

The solution to this problem was first formulated by Balitsky fl7 |. The Operator 



Product Expansion of the correlation functions of the type (|31| ) has a twofold structure. 
Terms of the first type come from the region x 2 ~ y 2 ~ 1/— p 2 and are proportional 
to vacuum expectation values (VEV) of local gauge invariant operators, multiplied by 
coefficient functions depending on p 2 , p\. In the following we refer to these terms as to 
local power corrections (LPC). 

Terms of the second type are called bilocal power corrections (BPC) and correspond 
to the contributions of large y 2 ^> x 2 ~ l/p 2 . In order to treat these terms for arbitrary 
q in the "Bjorken limit" p 2 ~ q ■ p — > oo one should expand the T-product of nucleon 
currents 

T[r ] (x)r ] (0)] = Y,C* L (x 2 )S n (x;0), (38) 

n 

in a series of non-local, gauge- invariant "string" operators of increasing twist n S n (x; 0), 
cf. 0. This expansion can be inserted into the correlation function ([[!]), producing a 
power series in 1/p 2 with coefficients given by the correlation functions of two nonlocal 
light-cone operators 



J dyexp (iq ■ y)(Q | T 



o u ' d (y + ^;y-^)s n (x;0) 



Q). (39) 



In general, calculation of the correlator (^) requires construction of a specific sum rule 
and may be very complicated. Remarkably, we have found that the most important BPC 
of dimension 6 can be evaluated exactly (i.e. related to the quark condensate) by using the 
equations of motion. The derivation essentially uses the Ward identity which we obtain 
in the Appendix. 

It should be noted that the OPE for the correlation function ( j3ll) for q 2 = is given 
by the sum of both LPC's and BPC's. In general only this sum has a physical meaning 
and is regularisation scheme independent. 
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Apart from these general remarks, we shall not go into details of the calculation which 
is relatively straightforward. A few more remarks are necessary, however, concerning 
the specific techniques used in the QCD sum rules approach to suppress contributions of 
higher states and further taking them into account in the duality approximation. 

To this end, the Borel transformation is applied to both sides of the sum rule, improv- 
ing the convergence of the operator product expansion series and suppressing exponen- 
tially contributions of higher resonances. The rationale for keeping nonzero value of the 
momentum transfer q in the above discussion is that in this kinematics one can consider 
p\ = p 2 and p\ = (p + q) 2 as independent variables and perform the Borel transforma- 
tion in both momenta. The advantage of this procedure is that in the double dispersion 
relation 

P(si,s 2 ) 



IT = / / dsids 2 -, 2\7 2\ + • • • ' ( 40 ) 

Jo Jo (si -p{){s 2 -pi) 

where p{s\,s 2 ) denotes the spectral density, one can ignore subtraction terms and con- 
tributions corresponding to non-diagonal transitions with singularities in only one of the 
two variables. In the standard duality approximation higher resonances and the contin- 
uum contribution are taken into account by the following model for the physical spectral 
density: 

p( Sl , s 2 ) = X 2 N 5( Sl - M 2 N )5(s 2 - M 2 N ) (N, Pl | 6 u ' d | p 2 , N) + Q( Sl - s )Q(s 2 - s )p c ( Sl , s 2 ). 

(41) 

where p c (si,s 2 ) is the corresponding spectral density calculated in perturbative QCD. 
Thus, by assumption, subtraction of the continuum contribution corresponds to con- 
straining the integration region in (si, S2)-plane to the duality region s±, s 2 < s . 

In the theoretical part of the sum rule the double Borel transformation is performed 
using the following general formula 

\[-vp 2 -vp 2 Y\ \ M 2 + M 2 ) v ; 

where t denotes the symmetric combination of Borel parameters M\ and Mf 

M 2 M 2 

t = M 2 + Mf ( } 

with v — 1 — v. In the symmetric case M 2 = Mf = 2t the subtraction of the continuum 
contribution corresponds to the replacement 



t n -> t n E n (t,s )] 

EJt,s ) = ( 1 - e' So/t , 

V ' 0; V L 1 (n-l)\\t 

in all terms of the OPE containing positive powers of the Borel parameter t. 



1 + — + ... + 



(44) 
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An explicit calculation leads to the following sum rule for the valence quark distribu- 
tions: 

e ~ Myt ^{^j } = £&CM, so) + f: >d (z,t, s ){(a s MG 2 ) + f%> d (z,t){qq) 2 

+ fs' d (z } t)(qq)-(qagGq} + ... (45) 

where we keep contributions of perturbation theory and operators up to dimension 8. 
The coefficient functions are given by diagrams shown in Fig. |6|. For the bilocal power 
corrections, examplified in Fig. [j], we take into account all contributions related to contact 
terms, and neglect contributions of correlation functions involving explicitly the gluonic 
fields in addition to quark fields which arise from the first term in the Ward identity in 

We find: 

/o"0) = -^t 3 E 3 (t, s ) f 1 du(9uu 2 + u 3 ) cos (uz) 

/o 0) = t!—/ E ^ s °) t M3uu 2 + u 3 ) cos (uz) (46) 

l07T 4 JO 

If 1 1 

f"( z ) = — — tEi(t, s ) / du[U8(u) + 6u-l + -z 2 u 3 ) cos (uz) + zu 3 u^ sin (uz)} 
96n d Jo 3 

1 r l 1 

fi( z ) = otEiCtjSo) / du[(25(u) + 3u — u H — z 2 u 3 ) cos (uz) + zu^u" 1 sin (uz)] 

9o7r d Jo 3 

(47) 

fe(z) = \ 

fUz) = ^cosz (48) 

r^) = —t- 1 

2 19 

f*(z) = r 1 (--cosz + —zsmz) (49) 



To keep the correct normalization of parton densities Q3"7|) the coupling in (fK3) should 



be substituted by the corresponding sum rule [16| 



e " A/ ArA,X^ = —L-t 3 E 3 (t, s ) + -^^tE 1 (t, So)((a s /Tc)G 2 ) + '^(qq) 2 -'^t~ 1 (qq}-(qo-gGq) + . . . 

(50) 

keeping the same terms in the OPE and using the same values of the Borel parameter 
t and the continuum threshold So as in (|4"5"D. Note that the coefficients in front of the 
dimension 8 terms in (49) and ( |50|) differ slightly from the corresponding ones in |15|] and 



TBI . The reason is that the authors of [T5|, [TBI evaluate the vacuum expectation values of 
nonlocal operators such as ^> (0)^> (x)^> (0)^> (x) reducing them to ~ ((^ (0)^> (x))) 2 using 
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the hypothesis of the dominance of the vacuum intermediate state, while we use this 
hypothesis in the calculation of local operators, like D 2 ^^?^ only. This leads to extra 

imfl. These terms are, however, suppressed as 



terms as compared to [15 



Nc being the 

number of colors, and the difference is not important numerically 

Please note that all coefficient functions in the sum rule in ( |4"5] ) are smooth functions 
of the Ioffe time, in contrast to the sum rules for the momentum fraction distributions 
given in Ref. |TJJ which contain expansions involving singular functions. Indeed, making 
a Fourier transform of our expressions, we obtain the following sum rule for momentum 
fraction distributions of valence quarks: 

cu.d 



-Mj,/t)2 _u,d, 



^nQv 



u 



fpert{u,t,s ) + f4 a {u,t,s )({a s /7c)G' 2 ) + fs' a {u,t){qq) 2 

+ is ,d K t)(gq) ■ (qvgGq) + ■■■ 



with the coefficients: 



f d M 

fl{u) 



(51) 



fo 


:«) = 


^t 3 E 3 (t, s )(9uu 2 + u 3 ) 


fo 


» = 


-^-t 3 E 3 (t, s )(3uu 2 + u 3 ) 

107T 4 


ft 


:«) = 


-^ 1 tE 1 {t,s )[{35{u)+6u- 


ft 


» = 


967T d 



(52) 



£*<«) 

2 1 Q 

-t-\-5(u) + -5'(u)) 



(53) 



(54) 



(55) 



where for any test function f(u) 



du 



tr 



/(„)= /'d^ (/(„)_ /(o)-u/ ; (0)) 



(56) 



It is easy to see that in high orders of the OPE the series of power corrections to the sum 
rules in coordinate and in momentum space will have the following typical behaviour 



t n n\ 



(AT 
t n n\ 



u 



t n n\ v 1 t n n\ v 1 



(57) 
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where (A) is the typical scale of vacuum fluctuations, of order (several hundred MeV) 2 , 
and the n\ suppression is due to the Borel transformation. For this reason the sum rules 
for Ioffe-time distributions can be justified theoretically at small values of z, while in the 
momentum space one faces a problem of the summation of the series containing singular 
functions. 

Let us now elaborate on this point. The calculation of parton distributions in coor- 
dinate space is essentially on the same theoretical footing as the calculation of moments 
as matrix elements of local operators. Mathematically, information about the moments 
is coded in derivatives of the coordinate space distributions at z — 0. Provided OPE con- 
verges fast enough, the ^-distribution is well defined and can be calculated at sufficiently 
small z by present nonperturbative approaches to QCD. In the particular technique of 
QCD sum rules, the results of calculations are usually assumed to be reliable provided 
contributions of vacuum condensates are sufficiently small, say stay within 30-40% of the 
total. It turns out that for valence quark distributions this criterium is satisfied for z < 3. 
We shall see that in practice the sum rule for u-quark distributions works in a larger 
interval, and for d-quark distributions for a shorter interval; for other approaches — e.g. 
lattice calculations — the limitations can be different. Our point is that for sufficiently 
small z one does not need to invoke any additional assumptions. The contributions of 
dimension 6 in fl49|) are respectable smooth functions at small z and must be taken into 
account, independent of their bad behavior at large z which produces 5- functions after 
the Fourier transform. 

In momentum space, the calculation of parton distributions pointlike in the Bjorken 
variable applies much more severe requirements to non-perturbative techniques, and in 
practice requires additional assumptions. In particular, the approach of Ref. |l5j assumes 
that singular terms in the OPE do not affect calculation of parton distributions at in- 
termediate values of momentum fraction u, and thus in this region all singular terms in 
the OPE can be neglected altogether (see also This would be true if summation of 

singular contributions produces a rather narrow smooth function with the support either 
in u ~ or in u ~ 1 regions. Our task in this paper is not to critisize this particular 
assumption, but rather to make clear that assumptions of this kind are always necessary 
to deal with parton distributions in momentum space, and thus provide an additional 
input. 

In fact, the assumption of Ref. |L5| is non-trivial, and to our opinion requires a 



better justification than given there. Mathematically, the statement about calculability 
of coordinate-space distributions at sufficiently small z does not imply calculability of 
momentum-space distributions at intermediate values of u. We find the neglect of singular 
contributions disturbing, since they are 100% essential for calculation of the moments, 
see Refs.fll^, [L4J. Since the calculation of the parton distributions along the lines of Ref. 
|TH| (and this paper) is only justified as the analytic continuation from the corresponding 
calculation of the moments 0, it is difficult for us to imagine that important contributions 
to the moments of the structure function will not show up in the distribution itself. 
Physically, the assumption about the small smearing of singular contributions implies 
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existence of a certain second scale in the hadrons, affecting the momentum distributions. 
We feel that a further study of this question is necessary, to prove that smearing of singular 
contributions does not affect the whole region of the Bjorken variable. Again, we repeat 
that advantage of coordinate space formulation is that it avoids making any assumptions 
of this kind, since singular contributions to not appear. 

Let us proceed to the description of our results. In the numerical analysis we use 
standard values of the parameters accepted in the QCD sum rules for the nucleon i.e., 
t ~ 1 GeV 2 and s = (1.5 GeV) 2 , and the following values for the condensates (the 
normalization point 1 GeV is implied): 

(qq) = (250 MeV) 3 , 
((a s /n)G 2 ) = 0.012 GeV 4 , 

(qagGq) = m 2 (qq); m 2 = 0.64 GeV 2 , (58) 

which correspond to the standard ITEP values rescaled to the normalization point /i 2 ~ 
M 2 N ~ 1 GeV 2 . 

The QCD sum rule prediction for the valence u-quark Ioffe time distribution in the 
proton U vaj \(z,fi 2 ~ 1 GeV 2 ) is shown as the thick solid line in Fig. |] and compared 
with an "experimental" distribution. The latter has been obtained from the leading-order 
parametrisation of Gluck, Reya and Vogt [TJ| normalized at 0.5 GeV 2 < /i 2 < 1 GeV 2 . We 



find remarkable agreement up to rather large values z < 4, corresponding to longitudinal 
distances in the proton rest frame of nearly 2 fm! When this QCD sum rule result is 
augmented by the assumption that U va \(z) is a sufficiently smooth function, and combined 
with large-z behaviour implied by the Regge theory ( see Eq.(Q)) it allows for complete 
reconstruction of U va \(z) and therefore of the distribution function. Lines marked as (a), 
(b) and (c) illustrate the relative importance of different contributions to the sum rule, 
and are obtained keeping in (|45|), ( |50"D the perturbative terms only (a), adding the gluon 
condensate contribution (b) and adding in addition also the (qq) 2 terms (c). Note that 
the U va x(z) distribution decreases at large z more slowly than the perturbative prediction, 
which is mainly due to bilocal corrections arising from large distances in the t-channel. 
The latter can be calculated as contact terms. We remind that these terms are discarded 
altogether in the approach of Ref. ]15| . 



Figure £5] shows stability of our prediction when the Borel parameter t is varied between 
1 GeV 2 (upper curve) and 1.5 GeV 2 (lower curve). One may conclude that in the region 
z < 4 the valence u-quark sum rule converges very fast and it is numerically stable. 

A sum rule similar to the one given by Eq. (|45| ) can be written for a non-zero value 
of the momentum transfer q 2 in the t-channel, allowing to study the radius of the valence 
quark distributions (cf. |2{|). We have checked that the radius of the valence u-quark 
distribution obtained from the sum rule in ([45]) is close to the measured electromagnetic 
radius of the proton, which is encouraging. 

The situation is not so good, unfortunately, for the valence d-quark distribution, see 
Fig. |It]. The sum rule prediction for D va [(z), shown as the thick solid line, is rather 
far from the leading-order GRV |18|] parametrisation, and is much more short-range. As 
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in Fig. |8|, we also show contributions of various terms in the OPE to the final result. 
It is seen that taking into account the gluon condensate contribution (b) improves the 
prediction considerably compared to the perturbative result (a), but this tendency is 
destroyed by the contributions of dimension 6 (c) and dimension 8. In the case of the 
d-quark distribution the bilocal corrections corresponding to the contact terms are absent, 
and the problem arises because of the diagram in which all nucleon momentum is carried 
by a single quark, see Fig. |6] (d), which contributes a term ~ cos(z) flUp. In the momentum 
space this contribution is proportional to 5(1 — u), where u is a fraction of longitudinal 
momentum carried by the quark, see fl55|) . Another reason is the absence of the bilocal 
power correcton of dimension 6, contributing a term ~ 5(u) in (|55|), suggesting that terms 
of higher dimension ~ S(u), 5'(u) . . . can be important. 

It turns out that the d-quark distribution is very sensitive to the numerical values of 
condensates of dimension 6 and 8 because of strong cancellation between corresponding 
terms. If we took qq = - (240 MeV) 3 and wig = 0.8 GeV 2 , which correspond to the normal- 
ization point fi 2 = 0.5 GeV 2 , we would obtain perfect agreement with the experimental 
analysis up to values of z of order of 1, where the sum rules prediction would abruptly 
turn down. Such a strong normalization-point dependendence makes our prediction for d 
quarks less reliable. 

A favourable structure of the OPE for the u-quark distribution and the complications 
for d-quarks have, presumably, no physical relevance, and are due to the particular struc- 
ture of the interpolating current fl33|). This choice is standard, but, as it follows from our 
analysis, not very convenient for the study of the d-quark distributions, since it implies 
that the correct behaviour of this distribution is due to higher order power corrections. 

A more detailed analysis of this problem goes beyond the scope of this paper. One 
could try a different interpolating current to improve the results for the d-quark distribu- 
tions, or calculate radiative corrections to the sum rule, which generally tend to soften the 
parton distributions (i.e., make them more extended in the Ioffe time) and are expected 



to be especially important for d-quarks, see [15 



However, it is worthwhile to demonstrate, at least semiquantitatively, that the higher 
order contributions indeed tend to smoothen the ~ 5(1 — u) contribution of dimension 6 in 
the d-quark sum rule, and are potentially able to bring it to the agreement with the data. 
To this end we use the concept of non-local condensates, introduced in jn| , which allows 
to consider the effects of the final correlation length in the QCD vacuum, the property 
that is missing in the local operator product expansion. 

Note that the contribution of the diagram in Fig. (d) is essentially proportional to 
the vacuum expectation value of the nonlocal operator (u aT (x)C^^u b (x))(u a (0)C'j^u bT '(0)) 
which produces the expansion 

((u aT (x)C^u b (x))(u a (0)C^u bT (0))) = ~(qq) 2 [l + ^m 2 x 2 + . . .] (59) 

where we have assumed the factorization to evaluate the coefficients. It is the expansion 
into the sum of local operators that generates the series of power corrections proportional 
to derivatives of 5(1 — u) in the sum rule for the momentum fraction distributions. As 



16 



noted in pifl , this expansion misses an important property of the correlation functions 
in Euclidian space, which is in existence of the final correlation length in the physical 
vacuum. To illustrate this point, let us consider the exponential parametrization 

9 roo 2 

((u aT (x)C^u b (x))(u a (0)C^u bT (0))) = -^(qq} 2 J o dve~ v f{v) (60) 

Moments of the function f(v) are determined by vacuum expectation values of local 
operators. The first few of them in the factorization approximation are fixed to be 



/ dvf{v) = 1 

JO 



duuf{u) = 2^T=2 m °' (61) 

On the other hand, one generally expects that the correlation functions in QCD decrease 
exponentially in Euclidian space, suggesting that 

((u aT (x)C^u b (x))(u a {0)C^u w (0))) ~ expf-Mu • v 7 ^ 5 ] (62) 

at large x 2 — > — oo, where Mr, is the correlation length which, loosely speaking, may be 
associated with the diquark mass. It is easy to see that this behaviour corresponds to the 
asymptotics 

f(y) ~ e ~ M > (63) 

at small v. Note that the expansion into the sum of local operators corresponds to the 
expansion of j[y) in delta-functions at v = 0. 

The effect of using the nonlocal condensate in the sum rules is easy to evaluate. In 
coordinate space, insertion of (|60D amounts to the substitution of the coefficient function 
fl in © by 

2 r 1 

fNLc( z ^) = o(<?<?) 2 ^ / duuf(ut) cos(uz). (64) 



3 

In momentum space this replacement is simply 



#(«) = !*(«) - f d NLC (u,t) = ^uf(ut) (65) 
For numerical estimates we choose a simple model 

m = y^^' m2Ju (66) 

with two parameters Mp and a. Equation (BT1) leads to the constraint 

M 2 

a - 3 = 2^ (67) 
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so that we only need to specify the correlation length. Its value has a direct physical 
meaning and is related to the difference between masses of heavy baryons, containing b 
quark and the uu pair, and the mass the h quark.| 



In Fig. 11 we show the r.h.s. of (pq) as a function of the momentum fraction it for 



the Borel parameter t — 1 GeV 2 and for two choices of Mr, = 700 MeV and 1 GeV. We 
see that the main effect of the finite correlation length is to push nonperturbative effects 
~ (qq) 2 away from the region it ~ 1. This may seem contrary to the physical intuition, but 
not necessarily so, because the effects of the quark condensate are qualitatively similar 
to the introduction of the constituent mass. We remark that this picture contradicts 
the expectations of Ref. f[7) |: Summation of singular contributions not only does not 



produce a narrow function with a support concentrated at u — > 1, but, on the contrary, 
the nonperturbative contributions die away at u — > 1 faster than any power of 1 — it. 
We stress that this is a direct consequence of of the final correlation length in the QCD 
vacuum. Our model estimates presented in Fig. [11] show that the resulting contributions 
are important at least up to u ~ 0.6. 

Numerical results for the valence d quark distribution D v (z) with the nonlocal con- 



densate are shown in Fig. |12[ We see that the situation improved somewhat at z > 2, 
although this type of contributions alone is not able to restore the agreement to the data. 
A better agreement can be obtained by choosing a larger value of the mixed condensate 
parameter wiq, but this possibility is not very attractive. We expect, however, that the 
results will be substantially improved by taking into account radiative corrections. An 
inspection shows that the sum rule for d quark distributions can also be saved by a large 
bilocal power correction, contributing a term ~ z 2 in coordinate space (alias ~ 5"(u) in 
momentum space). These corrections are difficult to evaluate, however, and we do not 
attempt this task in the present paper. An experience of QCD sum rule calculations gen- 
erally suggests that if there are indications that the sum rule is affected by contributions 
of high order in the OPE, it is advisable to use different interpolating currents for the 
participating hadrons. 



4 Summary and Conclusions 

We suggest to use Ioffe time distributions - the distributions of invariant longitudinal 
distances z essential in a deep-inelastic scattering process - as a suitable alternative to 
the conventional description in terms of momentum fraction parton distributions. The 
advantage of this formulation is that contributions of large and small longitudinal dis- 
tances that correspond in fact to different physics, become in this approach separated. 
The large-distances of distributions in coordinate space are governed by the Regge theory, 

§ We do not discuss this issue in detail, referring to a well-known relation between the asymptotics 
of the quark propagator at large distances in Euclidian space, and the difference between heavy meson 
and quark masses in the heavy quark limit |2^] . This difference is usually denoted by A and is one of 
the main observables in the heavy quark effective theory. For baryons the situation is quite similar. The 
range of values of Md used here corresponds to the estimates found in the literature. 
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and can be taken as an input, while the calculation at moderate distances can be within 
reach of current nonperturbative approaches, e.g. the lattice QCD. We illustrate these 
advantages using the QCD sum rule technique to calculate valence quark distributions, 
and compare it to the corresponding calculations in the momentum space. Although the 
sum rules derived in this paper are rather preliminary and can be improved significantly 
by taking into account further corrections, we obtain a very good description of the u 
quark distributions. The results for the d quark distributions are worse, the reason is an 
unfavourable structure of the operator product expansion series in this case, where higher 
order terms are important. 

We believe that our results can be improved by making a state-of-the-art QCD sum rule 
analysis, or with different techniques. In particular, we expect that Ioffe time distributions 
may be feasible for lattice calculations and for instanton models of the QCD vacuum of 



the type suggested in Ref. [23]. It would be most interesting to constrain in this way the 



polarized gluon distribution which is poorly known. 
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Appendix A 

In this Appendix we derive a Ward identity for the correlation function (|31~D useful in 
practical calculations. In addition, it makes normalization properties of valence quark 
distributions explicit. 

In our special kinematics vectors n M and are proportional: 

= ^—^ q (A.l) 
q ■ x 



Replacing 7 M n^ in the definition (|3^) of O l by ^ JnQ^ and integrating by parts over d A y 
in pi]) we get: 

n ■ x 



IT = — i j d A x exp (ip ■ x) J d 4 y exp (iq ■ y)- 



q ■ x 

9 -(Q | T{ v (x)¥(y + %,[y + %y- £]**(y - ^(0)} I 0) 



+ (A -> -A) . (A.2) 
The next step is to evaluate the derivative d/dy^ explicitly. One gets 

9 Wis A n r A A, Ti . A. 

-V\y + -)7 M [y + -; y - - -) = 



d Vtl y » 2"^ 2'" 2 J y » 2 
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D¥{y + |)[y + |; y - |] **(y - |) + f % + |)[y + |; y - j]L>¥(y - j) 

-f £ dmv + j)[y + f ; v + efb^G^y + e|)[y + 4^ - - y)> 

(A.3) 

where G a p is the gluon field strength. 

Inserting the right-hand-side of (|A.3|) into (|A.2| ) we note that terms containing D^f 



and lead to contact terms. To this end it is convenient to have in mind the functional 
integral representation of the correlation function in (|A.2| ), and use the identity 

exp (i j C d 4 x^j m(y) = ~^^y ex P (* J £ ^) ( A - 4 ) 

where C is the QCD Lagrangian. Making an integration by parts in the functional integral 
we obtain the delta function 5(x — y + ^) which allows to perform the integration over 
variable y explicitly. The net result can be written as a Ward identity: 

\ f f Tl ' X ( — ■ A 

IT = - / d 4 x exp (ip ■ x) dy exp [iq ■ y) / d£ (ft | Tli](x)^ l (y H ) 

2 1/ J c£ x J — 1 2 

x [y + f ; y + ^}gi a ^G aP { y + ^)[ y + tj-,y- - f W)} I ft) 



+ jV*z / (fx exp (ip • x) exp (ig • x)(ft I T(rf (x; - A)^(O)) | ft) 

j g • x L 

- (ft I T(7/(x)7f (0; A)) I ft)] + (A -> -A) (A.5) 

where N l = 1 and 2 for d- and u-quarks, respectively. The nonlocal currents r] l (x; — A) 
are defined as 

7] d {x; -A) = e abc u a (x) T C^u b (x)j 5 j"[x; x - A] c/ d / (a; - A) 

r] u {x; -A) = e abc u fT {x - A)C 7m [x - A; x] / V(x) 75 7 At cf (x) (A.6) 

The outcome of Eq. ( |A.5[ ) is that the complicated correlation function in ( [5T] ) is written as 
a sum of several simpler terms. The last two terms contain one integration less compared 
to the original expression, and therefore they cannot generate bilocal power corrections. 
The first term contains explicitly a gluon field. Thus, the OPE for this term starts with 
higher orders in the coupling and (or) the dimension of the corresponding operators. 
By an explicit comparison of the OPE applied to the correlation function (|3T| ) and to its 
equivalent form in ( |A.5| ) one can make sure that several important BPC's are transformed 
in this way to the LPC's related to vacuum expectation values of local operators, see 
Fig. [13] for the illustration. In fact, the last two terms in ( |A.5|) presumably collect all 



BPC's which can be reduced to vacuum condensates by low-energy theorems. Experience 
of practical calculations in the sum rules shows that these terms are typically the most 
important ones numerically. In this paper we do not take into account additional BPC's 
from the expansion of the term with an extra gluon in (|A.5|) . 
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In the limit A — > the last term disappears and we are left with the correlation 
function 

f 71 • T 

IT = NH / d 4 xexp (ip ■ x) (exp (iq ■ x) - 1) (ft | T (r](x)fj(0)) | ft) . (A.7) 

J q ■ x 

Equation (|A.7| ) can further be simplified writing 

(exp (iq ■ x) — 1) = iq ■ x / dv exp (ivq ■ x) , (A. 8) 

Jo 

so that we finally obtain 

IT = i\T du rv^-z / d A x exp (zp • a;) (ft | T (77(^)77(0)) | ft) 

= N l f'dvn^U^ip), (A.9) 
Jo ^dpf, 

where p = p + vq, arriving at the equation in fl35|). 

As an illustration of the use of the identity in (|A.5|) in practical calculation, let us 
consider the leading perturbative contribution (PT) to the correlation function fl34|) for 
u quarks, shown in Fig. §(a). First we note that owing to the explicit presence of the 
coupling in the gluon field, the first term in ( |A.5| ) can to this accuracy be neglected. Since 

(ft | T( v u (x; -A)tj(O)) I Q)pt = (ft | T( V (x)fj u (0; A)) | Q) PT , (A.10) 

the remaining contribution can be rewritten as: 

C 71 ' T 

W PT = N u i d 4 x exp (ip ■ x) (exp (iq ■ a;) - 1) (ft | T(r] u (x; -A)r/(0)) | ft) PT • 

(A.ll) 

and further using ( |A.7| ), as 

U u PT = N u iJ^dv J d 4 x exp (ip ■ x)(n ■ x)(Q \ T(if(x\ -A)r/(0)) | Vt) PT . (A.12) 

A straightforward calculation of (|A.12|) using dimensional regularisation produces 
1 1 r 1 T(3 — d) r 1 

-TriOl u = — 7 (p-n) 2 dv- ± 5^-7/ du (9uu 2 + u 3 ) cosuz , (A.13) 

4 4tt 4V ' Jo [-vp\-vplf- d Jo 

Comparing with equation ([34]) and performing the Borel transformation according to (|4 
we finally arrive at the u-quark coefficient function in equation ([47|). 
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Figure captions 



Fig. 1 Ioffe time distributions for partons at the scale /i 2 = 4 GeV 2 . The solid line denotes 
Gliick, Reya, Vogt (GRV) || set of parametrizations, the dashed line denotes CTEQ 
|| parametrization. U va \(z) and D V3 i(z) are the valence quarks distributions, Qu{z) 
and Qd{z) are the C-even up and down quark distributions, and G(z) is the gluon 
distribution. 

Fig. 2 Model of the Ioffe time valence-type distribution corresponding to the simple 
ansatz qv{u) = Nu a (l — u) 13 with a = —0.5 and (5 = 3. The normalization is such 
that Jq 1 duqviu) = 1. The dashed line shows the asymptotic expansion (0). Note 
that it matches almost perfectly the true behaviour for z > 6. 

Fig. 3 Model of the Ioffe time distributions for polarized gluon density Ag(u) = NGU a (l— 
u) 13 with a = and (3 — 4. The normalization constant Nq is chosen in such a way 
that the gluon polarisation Ag = 0.5. The two short-dashed curves were obtained 
by taking (3 equal to 3.5 and 4.5 respectively, and keeping a = and Ag = 0.5 fixed. 
The long-dashed curve is the asymptotic expansion (|2|). 

Fig. 4 Evolution of valence U v3b y(z,Q 2 ) and D vg ^(z, Q 2 ), and gluon G(z,Q 2 ) Ioffe time 
distributions. GRV f| parametrisation has been used. The solid line corresponds 
to Q 2 = 4 GeV 2 , the dashed line to Q 2 = 20 GeV 2 , respectively. 

Fig. 5 Evolution of up and down quark C-even Qu(z,Q 2 ) and Qd(z,Q 2 ) Ioffe time dis- 
tributions. GRV || parametrisation has been used. The solid line corresponds to 
Q 2 = 4 GeV 2 , the dashed line to Q 2 = 20 GeV 2 , respectively. 

Fig. 6 Typical diagrams contributing to the OPE of the correlator ([H]). 

Fig. 7. Generic form of a bilocal correction in the OPE of the correlator (|3"T|). 

Fig. 8 QCD sum rules calculation of the valence u-quark Ioffe time distribution function 
U va i(z, fi 2 ) at /i 2 ~ 1 GeV 2 . The thick solid line results from the OPE with operators 
with dimension 0, 4, 6 and 8. Thick dashed lines correspond to the leading-order 
QCD analysis of Ref.[18[] at the scales fi 2 = 0.5 and 1 GeV 2 , respectively. Solid 



line marked (a) is the perturbative contribution to the sum rule. Lines (b) and (c) 
describe respectively the sum rules with operators of dimension 4 and 6 taken into 
account. Note that VEV of dimension 8 gives already a small contribution. 

Fig. 9 Stability of the sum rule for U va \(z) against variation of the Borel parameter t. 
The upper curve corresponds to t = 1 GeV 2 , the lower to t = 1.5 GeV 2 . 

Fig. 10 QCD sum rules calculation of the valence d-quark Ioffe time distribution function 
Dvai{z, fi 2 ) at fi 2 ~ 1 GeV 2 . The thick solid line results from the OPE with operators 
of dimension 0, 4, 6 and 8. Thick dashed lines correspond to the leading-order QCD 



analysis of Ref.[18| at the scales jjtr = 0.5 and 1 GeV , respectively. Solid line 
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marked (a) is the perturbative contribution to the sum rule. Lines (b) and (c) 
describe respectively the sum rules with condensates of dimension 4 and 6 taken 
into account. 



Fig. 11 The nonlocal condensate contribution to the d-quark sum rule, see Eq.flHq). Lines 
marked as (a) and (b) correspond to the diquark mass parameters Mr, = 0.7 GeV 
and 1.0 GeV, respectively. Borel parameter t = 1 GeV 2 . 

Fig. 12 QCD sum rules calculation of the valence d-quark Ioffe time distribution func- 
tion D val (z) (solid lines). VEV of dimension 6 and 8 has been replaced by the 



phenomenological model of non-local four-quark condensate, Eq . (|60|) . Labels (a) 
and (b) refer to the diquark mass parameters M D = 0.7 GeV and 1.0 GeV, respec- 
tively. Thick dashed lines correspond to the leading-order QCD analysis of Ref 



normalized at /i 2 = 0.5 and 1 GeV 2 , respectively. The dashed line shows results of 
the standard sum rule (|45|). 



Fig. 13 Graphical illustration of the Ward identity, Eq. ( |A.5| ). Dashed lines denote path- 
ordered exponentials. 
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Figure 1: Ioffe time distributions for partons at the scale /i 2 = 4 GeV 2 . The solid line 
denotes Gliick, Reya, Vogt (GRV) || set of parametrizations, the dashed line denotes 
CTEQ ||] parametrization. U va \(z) and D va \(z) are the valence quarks distributions, 
Qu{z) and Qd{z) are the C-even up and down quark distributions, and G(z) is the gluon 
distribution. 
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Figure 2: Model of the Ioffe time valence-type distribution corresponding to the simple 
ansatz qv(u) = Nu a (l — u) 13 with a = —0.5 and (3 — 3. The normalization is such 
that Jo duqv(u) = 1. The dashed line shows the asymptotic expansion (0). Note that it 
matches almost perfectly the true behaviour for z > 6. 
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Figure 3: Model of the Ioffe time distributions for polarized gluon density Ag(u) = 
Ncu a (l — u) 13 with a = and (3 = 4. The normalization constant Nq is chosen in such 
a way that the gluon polarisation Ag = 0.5. The two short-dashed curves were obtained 
by taking (3 equal to 3.5 and 4.5 respectively, and keeping a = and Ag = 0.5 fixed. The 
long-dashed curve is the asymptotic expansion (|2J). 
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Figure 4: Evolution of valence U Y3 x(z, Q 2 ) and D v3u \(z, Q 2 ), and gluon G(z, Q 2 ) Ioffe time 
distributions. GRV || parametrisation has been used. The solid line corresponds to 
Q 2 = 4 GeV 2 ; the dashed line to Q 2 = 20 GeV 2 , respectively. 
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Figure 5: Evolution of up and down quark C-even Qu{z,Q 2 ) and Qd(z,Q 2 ) Ioffe time 
distributions. GRV || parametrisation has been used. The solid line corresponds to 
Q 2 = 4 GeV 2 , the dashed line to Q 2 = 20 GeV 2 , respectively. 
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Figure 8: QCD sum rules calculation of the valence u-quark Ioffe time distribution 
function U va \(z,fi 2 ) at fi 2 ~ 1 GeV 2 . The thick solid line results from the OPE with 
operators with dimension 0, 4, 6 and 8. Thick dashed lines correspond to the leading- 
order QCD analysis of Ref.|l8j at the scales /i 2 = 0.5 and 1 GeV 2 , respectively. Solid line 
marked (a) is the perturbative contribution to the sum rule. Lines (b) and (c) describe 
respectively the sum rules with operators of dimension 4 and 6 taken into account. Note 
that VEV of dimension 8 gives already a small contribution. 
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Figure 9: Stability of the sum rule for U va i(z) against variation of the Borel parameter 
t. The upper curve corresponds to t = 1 GeV 2 , the lower to t = 1.5 GeV 2 . 
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Figure 10: QCD sum rules calculation of the valence d-quark Ioffe time distribution 
function D va \(z,/j, 2 ) at fi 2 ~ 1 GeV 2 . The thick solid line results from the OPE with 
operators of dimension 0, 4, 6 and 8. Thick dashed lines correspond to the leading-order 
QCD analysis of Ref.fl8|| at the scales fi 2 = 0.5 and 1 GeV 2 , respectively. Solid line 
marked (a) is the perturbative contribution to the sum rule. Lines (b) and (c) describe 
respectively the sum rules with condensates of dimension 4 and 6 taken into account. 
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Figure 11: The nonlocal condensate contribution to the d-quark sum rule, see Eq. (|65|) . 
Lines marked as (a) and (b) correspond to the diquark mass parameters Mo = 0.7 GeV 
and 1.0 GeV, respectively. Borel parameter t = 1 GeV 2 . 
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Figure 12: QCD sum rules calculation of the valence d-quark Ioffe time distribution 
function D va \(z) (solid lines). VEV of dimension 6 and 8 has been replaced by the phe- 
nomenological model of non-local four-quark condensate, Eq . (|60"D . Labels (a) and (b) refer 
to the di quark mass parameters Mr, = 0.7 GeV and 1.0 GeV, respectively. Thick dashed 
lines correspond to the leading-order QCD analysis of Ref. [18] normalized at /x 2 = 0.5 and 



1 GeV , respectively. The dashed line shows results of the standard sum rule 
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Figure 13: Graphical illustration of the Ward identity, Eq. (|A.5|) . Dashed lines denote 
path-ordered exponentials. 
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